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ABSTRACT 

We present Kepler satellite photometry of KIC 10661783, a short-period binary star system 
which shows total eclipses and multi-periodic S Scuti pulsations. A frequency analysis of the 
eclipse-subtracted light curve reveals at least 68 frequencies of which 55 or more can be at- 
tributed to pulsation modes. The main limitation on this analysis is the frequency resolution 
within the 27-day short-cadence light curve. Most of the variability signal lies in the frequency 
range 18 to 31 d , with amplitudes between 0.1 and 4 mmag. One harmonic term (2 x /) 
O^l . and a few combination frequencies (/, + fj) have been detected. From a plot of the residuals 

versus orbital phase we assign the pulsations to the primary star in the system. The pulsa- 
tions were removed from the short-cadence data and the light curve was modelled using the 
Wilson-Devinney code. We are unable to get a perfect fit due to the residual effects of pulsa- 
tions and also to the treatment of reflection and reprocessing in the light curve model. A model 
where the secondary star fills its Roche lobe is favoured, which means that KIC 10661783 can 
be classified as an oEA system. Further photometric and spectroscopic observations will al- 
low the masses and radii of the two stars to be measured to high precision and hundreds of 
5 Scuti pulsation frequencies to be resolved. This could lead to unique constraints on theo- 
retical models of S Scuti stars, if the evolutionary history of KIC 10661783 can be accounted 
for. 

Key words: stars: binaries: eclipsing — stars: fundamental parameters — stars: oscillations 
— stars: variables: delta Scuti 



1 INTRODUCTION 

Eclipsing binary star systems (EBs) are our primary source of mea- 
surements of the properties of stars (Poppe rl ll980l; Undersell 199 ll ; 
iTorres et al.ll2010h . They are vital tracers of the physical processes 
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which govern stellar structure and e volution dAndersen et alll99Ct 
IPols et al.|[l997l ; lYoung et alj[200lh . From time-series photometry 
and spectroscopy of an EB it is possible to measure the masses 
and radii of the t wo stars empirically, and to accuracies of bet - 
ter than 1% (e.g. I South worth et all l2005bl : Iciausen etai]|2008h . 
These numbers in turn give the surface gravity and mean density 
of the stars. If measurements of the effective temperatures of the 
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stars are available, their lumin osities, absolute magnitudes and thus 
distance can be obtained (e.g. Guinan et al.ll 19981 : iHensberge et al.l 
2000; ISouthworth etail l2005ah . The accurately known surface 
gravities are also very helpful in performing a detailed chemi- 
cal abundance analysi s of the stars jPavlovsk i & South worthll2009l ; 
IPavlovski etaD 12009). The large number of physical properties 
measurable for EBs makes them prime objects for constraining the 
predictions of theoretical stellar models. 

8 Scuti stars are very promising targets for the study of stellar 
structure an d evolution through asteroseismology (see review by 
iBregedlioOOh . They are objects of spectral types A2 to F5 located 
near the main sequence inside the classical instability strip. Many 
show a large number of simultaneously excited non-radial pulsa- 
tion modes driven by the opacity mechanism. Photometric studies 
have revealed the presence of low-de gree (I 3) and low-radial 
order pressure modes (p-modes) jDaszvnska-D aszkiewi cz et alj 
2002), whereas spectroscopic studies showed that the se stars also 
oscill ate with high-degree modes up to £ — 14 ( Mantegazza 
|2004) . The presence of two unconnected convection zones allows 
the study of the efficiency of convec tive flux transport processes 
dDaszvfiska-Daszkiewicz et alj[200 5i). The abundance of gravity 
modes (g-modes) or mixed modes provides the opportunity to study 
regions close to the stellar core. 

Although there are several extensive observational studi es of 
8 Scuti sta rs from both the ground (e.g. lBregeretal.ll2005h and 
space (e.g. Buz asi et al . 2005), the successful application of astero- 
seismology has been hampered by uncertainty in their fundamental 
stellar parameters, an unknown mode-selection mechanism, com- 
plicated frequency spectra due to the presence of mixed modes 
and rotation, and the lack or uncertainty of the mode identifica- 
tion. The CoRoT space mission revealed an immense increase in 
detec ted frequencies in a 8 Scuti star compared to ground-based 
data dPoretti etal . 2009) but its interpretation was largely assigned 
to stellar gran ulation rather than all frequen cies being due to oscil- 
lation modes (Kallinger & Matthews 2010). The study of a 8 Scuti 
star in a binary system provides the possibility to accurately mea- 
sure its fundamental stellar parameters, greatly helping mode iden- 
tifications and thus the application of asteroseismology. 

One of the main hurdles facing studies of both pulsating stars 
and eclipsing binaries is the need to obtain extensiv e high-quality 
photometry. The recently-launched Kepler satellite jBorucki et al.l 
l20ld) provides an overwhelming solution to this problem. Kepler 
is currently monitoring roughly 150 000 pre-selected stars with ex- 
tremely high photometric precision and a duty cycle close to 100%. 
Almost all of these stars are observed in long-cadence mode, where 
the datapoints are sums of 270 consecutive 6 s exposures, leading to 
a sampling rate of 29.4244 min. Up to 5 12 stars at any one time can 
be observed instead in short-cadence mode, where sets of nine dat- 
apoints are summed to provide a sampling rate of 58.84876 s. The 
high duty cycle means that analyses of stellar pulsations are un- 
affected by aliasing (also termed spectral leakage) problems, and 
the high photometric precision leads to light curves of unparalleled 
quality. 

KIC 10661783 was found to be an EB as a result of the 
ASAS variability s urve}^] of the Kepler field jPoimahskilll997l ; 
IPigulski et ai]l2009h . We selected KIC 10661783 as a good target 
for the detection of pulsations in the component stars of EBs, and 
Kepler observations at both long and short cadence were obtained 
in the framework of the Kepler Asteroseismic Science Consortium 
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Table 1. Basic observ able properties of KIC 10661783; 

Refere nces: (1) Tycholttog et al.Ul997h : (2) 2M ASSlCutri et alJl2003l); ( 3) 
KIClO lBrown et alj <2005n ; (4) AGK3 catalogue iBucciarelli et all 19921) . 



KIC 10661783 


Q2000 


1921 11.619 


1 


52000 


+47 58 42.95 


1 


Tycho identification 


TYC 3547-2135-1 


1 


2MASS identification 


119211161+4758430 


2 


B T 


9.844 ± 0.022 


1 


V T 


9.563 ± 0.021 


1 


^2MASS 


9.264 ± 0.024 


2 


-H2MASS 


9.175 ±0.026 


2 


-f^2MASS 


9.177 ±0.022 


2 


Kepler magnitude 


9.586 


3 


Spectral type 


A2 


4 



(KASC: lGUliland et alJ[2oTobl> . The resulting light curves clearly 
exhibit deep total eclipses as well as multi-periodic pulsations with 
frequencies representative of short-period 8 Scuti pulsations. 

Using the 2MASS magnitudes o f the system and the InfraRed 
FLux Method l lBlackwell et alj|l980h we find an effective temper- 
ature of T eS = 8000 ± 160 K. The light contribution of the sec- 
ondary star will act to lower the T e g found by the IRFM, s o the pri- 
mary star will be a few hundred degrees hotter than this (ISmallevI 
1993). Some basic observable properties of KIC 10661783 are col- 
lected in TableQ] We find that it is a semidetached binary whose 
evolutionary history will have been strongly affected by interac- 
tions and mass transfer between the stars. This will complicate the- 
oretical analyses of the frequency spectrum of the pulsating pri- 
mary star, and may in the future require a new theoretical treatment 
of such objects. 

The study of pulsating stars in EBs has the potential to be a 
remarkably exacting test of theoretical stellar models due to the 
large number of physical constraints which can be applied. In the 
case of 8 Scuti stars, the measurement of a large number of fre- 
quencies in a star of known mass and radius should enable identi- 
fication of the pulsation modes. This previously elusive goal is an 
important step towards obtaining a proper physical understanding 
of 8 Scuti pulsations. In addition, it holds the promise of being able 
to directly detect the spatial distribution of pulsations on a stellar 
surface via eclipse-mapping techniques, leading in turn to empir- 
ical mode identificatio ns. Although few pulsating EBs are known 
dRodrfguez e t al. 2004), Kepler will lead to a major increase in this 
number. An alternative approach which is possible for the brightest 
nearby stars is interfer ometric observations to pin down the radius 
of a pulsating star (e.g.)Monnier e t al.ll201ol) . 



2 DESCRIPTION OF THE OBSERVATIONS 

KIC 10661783 was observed in both long and short cadence by 
the Kepler satellite. Detailed descripti ons of the characteri stics of 
these observations can be found in Ijenkins et all l l2010allbl) and 
iGilliland et alj yOlOa). The short-cadence observations comprise 
39 810 datapoints obtained over 27 days in Quarter 2.3, and the 
long-cadence data encompass 476 datapoints from Quarter and 
1639 from Quarter 1. Following the advice in the release notes ac- 
companying the data, the uncorrected flux measurements were used 
in our analysis. The timestamps of the data are expressed as BID 
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on the UTC timescale, and we converted the flux measurements to 
a relative magnitude scale. 

In order to extend the time interval covered by the available 
photometric data we obtained a l ight curve of KIC 10661783 from 
the Su perWASP public archiveQ dPollacco et alj200d : lButters et al.l 
l20ld) . This contains 9098 datapoints, obtained between 2004 
May and 2008 Au gust, which have been detrended using the 
iTamuz et alj d2005h algorithm. They were converted to a relative 
magnitude scale, transformed to the BJD(UTC) timescale, and 
3.5<r-clipped to reject outliers. We did not use the ASAS data as 
only 66 high-quality brightness measurements are available from 
this survey. 



3 ORBITAL PERIOD AND PRELIMINARY ECLIPSE 
MODEL 

In order to measure the orbital ephemeris of KIC 10661783 we fit- 
ted an eclipse model to th e light c urves using the JKTEBOP cod^fl 
dSouthworth etal] l2004allbl 120071) . JKTEBOP represents the stars 
as biaxial spheroids (Etzel 1981). This simple model is extremely 
quick to evaluate and has very low numerical noise, so is well suited 
to analysing data which are both numerous and of high precision. 

The Kepler short- and long-cadence data and the SuperWASP 
data were initially fitted separately, after which it was possible to 
bring them to a common magnitude scale as well as assign realistic 
measurement errors. An analysis was th en executed on the com- 
bined data, and Monte Carlo simulations dSouthworth et al.ll2004d, 
2005b) were performed in order to estimate the errorbars. The dif- 
ference in the response functions of Kepler and SuperWASP is 
small enough to have no significant effect on the results. From this 
process we find an orbital ephemeris of: 

Mini = BJD(UTC)2455065. 77701(5) + 1.23136220(24) x E 

where the parenthesised quantities give the la uncertainty in the 
previous digit. 

For all subsequent analyses we used only the short-cadence 
Kepler data, as these have a high precision, a duty cycle of 89%, 
and a sampling rate which is much faster than the photometric vari- 
ations seen in the light curve. JKTEBOP was used to provide a fit to 
the data using the orbital ephemeris determined above. Whilst the 
primary component of KIC 10661783 is deformed bey ond the lim- 
its of applicability of the EBOP model ( Popper] |l98l), the system 
is within the range where J KTEBOP is still capable of p roviding a 
good fit to the light curve (Southworth & Clausen 2007). Once a 
good morphological model was obtained, it was subtracted from 
the observations, leaving behind a residual light curve containing 
primarily the S Scuti pulsations and noise. 



4 THE PULSATION CHARACTERISTICS OF 
KIC 10661783 

We performed a frequency analysis of the residual light curve in 
order to search for pulsation frequencies. This light curve (see 
Figure [2} clearly shows variations on a timescale of hours with 
a total peak-to-peak; amplitude of about 30mmag arising from 
multi-periodic pulsations of at least one of the components of the 

2 SuperWASP archive: jhttp : / /www . wasp . le . ac . uk/public/| 

3 JKTEBOP is written in FORTRAN77 and the source code is available at 

|http : / /www .astro . keele ■ ac ■ uk/~ jkt/| 



binary system. We searched for pulsation frequencies in an it- 
erative way by making use of Lomb-Scargle periodograms and 
multi-pe riodic non-linear least -squares fitting algorithms using Pe- 
RIOD04 fcenz & B reger 200% and our own codes. 

At each step of the iterative algorithm we selected the 
frequency peak with the highest amplitude and computed a 
multi-periodic least-squares fit to the data using all detected 
n frequencies simultaneously by applying the formalism Z + 
Ai sin(27r/it + <f>i), where Z denotes the zero point, Ai is 
the amplitude of the frequency fi, t is the time of each measure- 
ment, and <j)i the phase. The data were then pre-whitened with the 
derived fit for further analysis. This algorithm was repeated until no 
significant frequency peaks were detected. For the first 50 frequen- 
cies we included each frequency in the least-squares fit. For the 
other 1 8 frequencies we kept all previous frequencies fixed at their 
known value, to avoid problems due to strong correlations between 
parameters. To determine the significance of a frequency peak, we 
adopted the signal-to-noise ratio cri terion S/N > 4 in amplitude 
dKuschnig etalll 1 9971 : iBregeJI 1 9931) . The S/N of all detected fre- 
quencies is listed in Table [2] 

The periodogram of the light curve (Fig. [3] top panel) shows 
that the main variability of KIC 10661783 lies in the frequency re- 
gion 20-30 d" 1 (231-337 /zHz). Additional signal is visible at 50- 
60 d _1 and at low frequencies (<5 d ). The low-frequency signal 
is at least partially due to imperfect removal of the eclipses from 
the light curve (/2s = / or b and its harmonics) and due to long-term 
trends that are present in the Kepler data (/?). We also checked for 
significant frequencies up to the Nyquist frequency at 680d _1 , but 
could detect no signal above 60 d -1 . 

Table[2]lists all frequencies that were detected in the data with 
an amplitude above S/N ^ 4. We measured the noise in regions 
of three different sizes around each frequency. Due to the large 
number of frequencies the noise level is certainly overestimated, 
which can be seen when looking at the noise level above 70 d -1 
where no stellar signal is present. Nevertheless, we believe that it is 
a good choice to stop with the frequency analysis at this stage since 
many of the remaining peaks ar e closer than the theoreti cal fre- 
quency resolution of 0.05 d _1 dLoumos & Deeming! 1 19781) to the 
already detected frequencies. Moreover, even small amplitude or 
phase changes of pulsation modes show up as additional frequency 
peaks in the Fourier spectrum. As we will show below, amplitude 
variations seem to be present in this star, at least in part due to the 
eclipses. 

At least 68 significant peaks were found in the data, with 
frequencies between 0.01 and 59 d _1 and amplitudes from 4 to 
0.1 mmag. The frequency range which can be associated with non- 
radial p-mode pulsation is 14— 34 d -1 , where we have detected 55 
peaks. Some of the pulsation frequencies in this region are sep- 
arated by only the theoretical frequency resolution (0.05 d -1 ) or 
less. At this stage we cannot decide whether these close frequen- 
cies are due to amplitude variations or pulsation modes intrinsic 
to the star. A longer time-series is required to solve this question. 
Clos e frequencies have been reported before in s everal S Scuti stars 
(e.g. lBreger & Bischofll2002l : lBreger et al.ll2005l) and have been at- 
tributed to the presence of mixed modes and modes having dif- 
ferent order and different ^-values. Since KIC 10661783 is already 
evolved, mixed modes could be the cause for the observed high 
number of modes between 20 and 30 d _1 . 

We examined if the peaks between 50 and 60 d" 1 are due to 
combination frequencies (fi + fj) or harmonics (N x /) of the 
terms between 20 and 30 d _1 . In this region, we found three fre- 
quencies at 52.39, 55.18 and 58.76 d" 1 . The highest peak in this 
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Figure 1. The full Kepler short-cadence light curve of KIC 10661783. 



region at 58.76 cP 1 can be identified as the harmonic 2 x fz which 
implies that fz deviates from a pure sinusoid. Surprisingly, such a 
harmonic term is not found for either /i or /2. Several different 
possible combinations exist for the other two high-frequency peaks 
and we have listed the most probable combinations in Table|2] We 
cannot exclude that the frequencies between 50 and 60 d _1 are in- 
trinsic pulsation modes. 

The lowest-frequency signals at /V and /si are probably due 
to long-term trends that are present in the data, rather than high- 
order g-modes. The time base of the light curve is too short to in- 
terpret those further. In addition, low-frequency signal is detected 
at seven frequencies which appear to be harmonics of the orbital 
frequency to within the frequency resolution. This holds the po- 
tentially interesting outlook of modelling them as tidally excited 
modes, such as recently detect ed in CoRoT photom etry of the ec- 
centric B-type EB HD 174884 iMaceroni et al.l2009h and in Kepler 
data of the h ighly eccentric spectacular A-type binary HD 187091 
(Welsh et alj|201lh . In the latter object, all detected frequencies 
after orbit subtraction occur at multiples of the orbital frequency, 
some at rotationally split values according to the rotation of the 
primary. This is exactly the signature expected for tidally excited g~ 
mode s and this detection was inteipreted as such ( lWillems & Aertsl 
l2002h . In the case of KIC 10661783, it is less likely that we are 
dealing with tidally excited modes, given that the orbit is circular. 
Nevertheless, in such a close binary, one coul d still be dealing with 
tidal p erturbations of intrinsic oscillations dRevniers & Smeversl 
l2003aH bh. Before embarking on such an interpretation, we first need 
to determine much more accurate values of the frequencies in the 



residual light curve, as well as of the fundamental parameters of the 
binary components, after excluding that these orbital harmonics are 
due to imperfect removal of the eclipses from the light curve rather 
than intrinsic oscillations. 

In order to check if the detected frequencies are consistent 
with those predicted by theory we computed the oscillation spec- 
trum of a typical stellar m odel for a mass of 2 Mg with the code 
CLES ( iScuflaire et al. 2008), using the input physics as described in 
iDegroote et al.N2010l) . We used an initial hydrogen fraction of 0.70, 
a metallicity of Z — 0.02, a convective overshooting parameter of 
0.2 local pressure scale heights and an age of 8 x 10 8 yr, which 
corresponds to T c g ~ 8000 K. The oscillation spectrum for modes 
of degree £ = 0, 1, 2 was computed with the pulsation code MAD 
(Dupret et al. 2002), taking into account rotational splitting for an 
equatorial rotation velocity of 95kms _1 , which we deduce from 
the model stellar radius (2.3 Rq) and assuming the rotation period 
to be equal to the orbital period. It is justified to work with an evo- 
lutionary model ignoring rotation for this consistency check, given 
that the rotational frequency is less than 20% of the cricital fre- 
quency in the approximation of co-rotation with the orbit. We find 
hundreds of eigenfrequencies in the measured range between and 
70 d _1 , where the gr-modes occur below 16 d -1 and the p-modes 
above 13 d~ 1 . It is thus easy to explain the measured oscillation 
spectrum in terms of p- and g-modes in a qualitative sense without 
having to invoke mode degrees above 1 = 2. We therefore con- 
clude that the detected frequencies are consistent with theoretical 
predictions of low-degree modes. At present our frequency resolu- 
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Figure 2. Light curve of KIC 10661783 after removal of the eclipses. The top panel shows the complete time series of Q2.3. Beating is clearly visible and is 
due to the presence of multiple pulsation frequencies. The bottom panel shows a short section of the light curve compared to a 68-frequency fit to the data. 



tion is insufficient to perform mode identification and exploit the 

freq uencies seismically. 

Kallinger & Matthewsl l20ld) suggested that the dense fre- 
quency spectra in A-type stars deduced from CoRoT data could 
be due to granulation rather than individual pulsation modes. The 
primary of KIC 10661783 does indeed follow the scaling law pro- 
posed in that study, but we have no reason to doubt that the fre- 
quencies we report here are due to individual oscillation modes. 



4.1 The effects of binarity 

To rule out that the many detected frequencies are caused by chang- 
ing amplitudes during the eclipses we also analysed the light curve 
after clipping out the times of the primary and secondary eclipses. 
This revealed essentially the same frequencies as the light curve 
including eclipses; the amplitude spectrum of the residuals also 
closely resembled the one in Fig.[3] We therefore conclude that 
the large number of frequencies is not related to the presence of 
eclipses. 

Since KIC 10661783 is in a multiple system one can expect 
light-time-travel effects which affect the observed pulsation phases. 
From the derived orbital parameters of the system we estimated a 
maximum light time effect of approximately 2.3 s. We examined if 
such variations are present in the data by following the procedure 



described in lBregend2005l) . We assumed that all observed frequen- 
cies originate from the same star (which is not necessarily true) 
and split the light curve into ten separate datasets phased with the 
orbit. For each dataset we determined the maximum of the cross- 
correlation of the data with the fit using all 68 frequencies. If the 
light curve is shifted in time due to the light time effect the peak 
in the cross-correlation would be expected at the value of the shift. 
The top panel in Fig.|4] shows the derived shift plotted against the 
orbital phase. No significant variation is visible, which may be due 
to the limited frequency resolution or size of the dataset. But we see 
a significant increase of the height of the residuals (bottom panel) 
at the times of eclipses which shows that the fit is worse at the 
times of the eclipses. This is because we see a different amount 
of light from the pulsating star during the eclipses compared to at 
other times. From Fig.|4]we see that the residuals are the greatest 
during primary eclipse, which indicates that the primary star is the 
pulsating component in this binary system 

To further test what the contribution of each star to the ob- 
served frequency spectrum is, we examined the light curve at the 
times of the eclipses. During secondary eclipse, the primary com- 
pletely blocks the light of the secondary. Any observed oscillation 
frequencies therefore arise only from the S Scuti star primary. Dur- 
ing primary eclipse, a partial blocking of the light of the S Scuti 
star should have an effect on the observed amplitudes of its non- 
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Figure 3. Amplitude spectra of KIC 10661783 before (top panel) and after pre-whitening with 15 (middle panel) and 68 frequencies (bottom panel). There is 
clearly signal left after the removal of 68 frequencies but the limited frequency resolution prevents the reliable detection of further peaks. Note that the y-axis 
scale varies for the different panels. 



radial pulsation modes. This can in principle enable mode identi- 
fication by comparing in-eclipse and out-of-eclipse pu lsation am- 
plitudes through eclipse mapping (e.g. Reed et al. 2005). Since the 
eclipses are regular and relatively short events, the frequency spec- 
trum shows both a much higher noise level and strong aliasing with 
peaks separated by the orbital frequency (see Fig. [5}. One therefore 
cannot expect to detect as many frequencies as for the analysis of 
the complete light curve. 

The top diagram of Fig.|5]shows the amplitude spectra during 
the primary (top panel) and secondary (bottom panel) eclipses be- 
tween and 80 d _1 . In both cases, the main signal lies between 
20 and 30 d _1 and there are no significant differences visible in 
the data before pre-whitening. The comb-like structure is due to 
the strong aliasing from selecting only data during eclipses. The 
spacing between adjacent peaks is equal to the orbital frequency of 
about 0.81 d~\ We then pre-whitened both datasets with the same 
ten frequencies with the highest amplitude. Of these ten frequen- 
cies, eight have also been found in the analysis of the complete 
light curve. The slight differences in the frequency values between 
the datasets (see Tables [3] and |4]l is due to the fact that we also im- 
proved the frequency values during the least-squares fitting. The 



term / a can be associated with a slow zero-point drift in the light 
curve, similar to previously found ft. The correct value of the term 
/b = 15.43 d~ , which has been found in both eclipse datasets 
but not in the complete light curve solution, is uncertain due to 
the strong aliasing. The fact that this term is only seen during the 
eclipses could be due to the complex shape of the spectral window 
and the overlapping of several aliasing peaks. 

The amplitude spectra of the pre-whitened datasets is shown 
in the lower diagram of Fig.|5] It is evident that during primary 
and secondary eclipse, the frequency content of the residuals dif- 
fers. During primary eclipse, when the cooler star is in front of the 
8 Scuti star, there is a lot of signal left below 20 d _1 in contrast to 
the secondary eclipse, when only the 8 Scuti star is visible. This can 
be due to the fact that some of the pulsation modes of the 8 Scuti 
star change their visibility due to the partial eclipse. If we assume 
that the rotation axis of the S Scuti star is perpendicular to the line 
of sight, the eclipse configuration derived in Section|5] would lead 
to an increase in the visibility of axisymmetric pulsation modes 
during eclipses. During secondary eclipse, there are still many fre- 
quencies left in the region between 20 and 30 d~ , which could be 
identified with terms that were detected in the complete light curve 
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Table 2. Frequencies in the order of detection in the amplitude spectrum 
of KIC 10661783. The noise was computed in a range of 5 d —1 around 
each frequency. The formal uncertainty in the last digit as derived from 
the least-squares fit is indicated in brackets. Since the frequency was held 
constant for the final fit of all 68 frequencies, we do not list the formal uncer- 
tainties of the frequency. T he theoretical frequency resolution is 0.05 d _1 
iLoumos &~D eeming 193). 



Table 2 - continued 



Designation 


Frequency 


Amplitude 


Phase 


.3 /IN 




(A — 1 \ 


(mmag) 


(rad) 




ft 


28.135 


4.163 (8) 


2.114(2) 


98.7 


h 


22.341 


3.114(8) 


1.351 (3) 


64.4 


ft 


29.383 


2.396 (8) 


1.306(4) 


59.2 


ft 


25.902 


1.745 (8) 


2.906 (5) 


32.7 


h 


27.810 


1.716(8) 


1.899 (5) 


39.2 


h 


24.408 


1.366 (8) 


1.156(7) 


26.7 


ft 


0.0166 


1.157(8) 


1.108 (23) 


31.5 


ft 


29.759 


1.022 (8) 


2.901 (8) 


27.0 


ft 


23.462 


0.972 (8) 


1.465 (9) 


19.5 


fto 


28.622 


0.839 (8) 


1.143 (10) 


18.8 


fn 


21.078 


0.754 (8) 


2.591 (11) 


16.4 


fl2 


25.433 


0.656 (8) 


2.155 (13) 


12.6 


fl3 


25.802 


0.633 (8) 


1.083 (13) 


11.7 


fl4 


27.063 


0.712(8) 


0.683 (12) 


15.1 


/is 


23.016 


0.567 (8) 


0.814(15) 


10.8 


he 


25.192 


0.571 (8) 


0.166(15) 


10.9 


fn 


29.442 


0.552 (8) 


2.064 (15) 


14.0 


fts 


28.655 


0.456 (8) 


2.859 (19) 


10.8 


fl9 


26.840 


0.485 (8) 


0.803 (17) 


9.3 


f%0 


27.001 


0.402 (8) 


1.452 (21) 


8.7 


ftl 


24.618 


0.414 (8) 


1.845 (20) 


8.1 


fl2 


27.860 


0.401 (8) 


2.004(21) 


9.4 


ftz 


23.890 


0.348 (8) 


2.824 (25) 


6.7 


fti 


20.718 


0.337 (8) 


0.320 (24) 


7.5 


ft$ = forb 


0.809 


0.313 (8) 


2.301 (27) 


9.2 


he = 2/3 


58.764 


0.313 (8) 


1.404(27) 


18.5 


ftt 


22.569 


0.320 (8) 


2.997 (26) 


6.4 


fts 


24.448 


0.293 (8) 


2.744 (29) 


6.3 


/29 


21.662 


0.301 (8) 


1.029 (29) 


6.1 


/30 


21.839 


0.285 (8) 


0.475 (29) 


6.2 


/31 


27.318 


0.279 (8) 


0.064 (30) 


6.2 


f-',2 


20.607 


0.279 (8) 


1.752 (30) 


6.1 


/33 


28.948 


0.273 (8) 


2.421 (30) 


6.9 


/34 


21.035 


0.263 (8) 


2.831 (33) 


5.5 


/35 


26.464 


0.263 (8) 


2.385 (32) 


4.9 


f'-'M 


20.438 


0.273 (8) 


2.426 (32) 


5.7 


f'',7 


27.537 


0.287 (8) 


0.396 (29) 


6.0 


/: J >8 


26.800 


0.243 (8) 


1.844 (34) 


5.1 


/39 


30.176 


0.249 (8) 


1.175 (33) 


6.9 


/40 = 5f or b 


4.062 


0.244 (8) 


2.528 (34) 


8.6 


hi 


24.283 


0.244 (8) 


1.541 (34) 


5.8 


h'2 


22.073 


0.224 (8) 


2.788 (38) 


4.7 


Us 


31.211 


0.204 (8) 


1.384 (42) 


6.5 


j ! 44 


25.402 


0.215 (8) 


1.889 (40) 


4.3 


fiS 


25.850 


0.230 (8) 


0.718 (37) 


4.3 


he 


22.166 


0.217(8) 


2.302 (39) 


4.4 


fvi = Zforb 


1.600 


0.221 (8) 


2.026 (38) 


5.8 


/t8 


22.793 


0.220 (8) 


2.669 (38) 


4.1 


/t9 


20.248 


0.219(8) 


0.201 (37) 


4.5 


ho 


27.779 


0.232 (8) 


1.200(38) 


5.2 


hi 


28.231 


0.200 (8) 


2.330 (41) 


4.7 


hi 


24.214 


0.210 (8) 


1.183 (39) 


4.3 


hs 


27.690 


0.193 (8) 


3.132(41) 


4.3 


f&A 


18.868 


0.182 (8) 


2.810(50) 


4.2 


f55 


24.380 


0.208 (8) 


0.229 (45) 


4.2 


he = if orb 


2.434 


0.171 (8) 


2.435 (49) 


4.8 



Designation 


Frequency 


Amplitude 


Phase 


S/N 




(d 1 ) 


(mmag) 


(rad) 




/ 57 


28.788 


0.181 (8) 


1.411 (48) 


4.0 


J DO 


30.403 


0.158 (8) 


1.620 (53) 


4.2 


joy j 1 \ .1 14 


55.189 


0.154 (8) 


2.162 (54) 


4.8 


/60 = §forb 


4.868 


0.157 (8) 


2.269 (53) 


5.6 


hi 


0.092 


0.166 (8) 


0.197 (55) 


4.5 


fe2 


17.236 


0.155 (8) 


2.313 (54) 


4.2 


fe:i = h + fts 


52.395 


0.152 (8) 


2.590 (55) 


4.0 


fe4 


31.159 


0.152 (8) 


2.085 (71) 


4.0 


/65 = 9 f or b 


7.306 


0.129 (8) 


1.628 (64) 


5.6 


he 


14.965 


0.116(8) 


1.960 (72) 


4.2 


hi 


12.402 


0.102(8) 


0.208 (77) 


4.6 


h8 = ll/orf) 


8.935 


0.108 (8) 


1.355 (78) 


5.7 


residuals 




0.00116 
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Figure 4. Top panel: O— C shifts as a function of the orbital phase assuming 
that all 68 frequencies originate from the same star. Bottom panel: Residuals 
of the data minus the 68-frequency solution at different orbital phases. The 
residuals are significantly higher at the eclipse phases. 



Table 3. List of frequencies that were detected by only looking at the phases 
of the primary eclipse (when the light of the S Scuti star is partially blocked 
by the secondary component). The frequency designation of the previously 
found frequencies in Tablets used. The formal uncertainty in the last digit 
is indicated in brackets. 



Designation 


Frequency 


Amplitude 


Phase 




(d- 1 ) 


(mmag) 


(rad) 


ft 


28.133 


2.84 (5) 


3.41 (1) 


ft 


22.339 


3.61 (5) 


1.72(1) 


ft 


29.388 


2.66 (5) 


3.02 (2) 


h 


25.901 


2.77 (5) 


0.20 (2) 


h 


27.811 


2.58 (5) 


2.65 (2) 


h 


24.415 


1.14(5) 


5.64 (4) 


fl2 


25.447 


0.94 (5) 


4.44 (5) 


fl3 


25.815 


0.95 (5) 


2.37 (5) 


/a 


0.023 


0.55 (5) 


4.68 (14) 


/b 


15.394 


0.91 (5) 


0.35 (5) 


residuals 




0.00162 
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Figure 5. Amplitude spectra of the light curve during the times of the 
eclipses. In each diagram, the upper panel shows the amplitude spectrum 
during primary eclipse and the inverted lower panel shows the amplitude 
spectrum during the secondary eclipse. Top diagram: Amplitude spectra be- 
fore pre-whitening. Bottom diagram: Amplitude spectra after pre-whitening 
of the same ten frequencies for both datasets. The comb-like structure is due 
to aliasing. 

(see Table|4j. The frequencies listed in Table[4]originate solely from 
the 8 Scuti star and can thus be used as the basis for seismic stellar 
modelling. 



5 LIGHT CURVE ANALYSIS 

Whilst JKTEBOP provided an impeccable fit to the Kepler 
light curve, the model parameters are not expected to be re- 
liable for stars which are as distorted as the components of 
KIC 10661783. We have theref ore fitted the short-cadence data us- 
ing the Wilson-Devinnev code {Wilson & Devinne-dfilml : [Wilson 
ll979tlWilson & Van Hammell2004) . which implements Roche ge- 
ometry to accurately represent the surface figures of stars in binary 
systems. We used the 2004 version of the code (hereafter referred to 
as WD2004), with automatic iteration performed using the JKTWD 
wrapper written by JS. The primary and secondary components are 
referred to as star A and starB, respectively. 

In order to save computing time we combined the light curve 
into 230 phase bins using the orbital ephemeris obtained in Sect.[3] 
with a finer bin size during the eclipse phases. The orbital period 
and time of primary mid-eclipses were then fixed at 1.0 and 0.0, 
respectively. The stars were set to rotate synchronously with the 



Table 4. List of frequencies that were detected by only looking at the phases 
of secondary eclipse (when only the 8 Scuti star is visible). The frequency 
designation of the previously found frequencies in Table|2]is used. The for- 
mal uncertainty in the last digit is indicated in brackets. 



Designation 


frequency 


Amplitude 


Phase 




( a ) 


(mmag) 


[ma) 




28.134 


3.07 (3) 


0.21 (1) 


h 


22.338 


3.56 (3) 


2.31 (1) 


h 


29.383 


2.67 (3) 


2.62(1) 


h 


25.902 


2.02 (3) 


2.76 (2) 


ft 


27.810 


0.91 (3) 


4.84 (4) 


fe 


24.407 


2.34 (3) 


4.31 (1) 


ho or /ig 


28.659 


1.00(3) 


5.74 (4) 


/u 


21.020 


0.85 (3) 


3.43 (4) 


fl2 


25.428 


1.36 (3) 


3.34 (3) 


ha 


25.808 


1.00(3) 


1.61 (3) 


he = 2/3 


58.765 


0.33 (3) 


1.93 (10) 


/63 


55.134 


0.29 (3) 


0.92 (11) 


far 


52.463 


0.37 (3) 


0.54 (9) 


/a 


0.026 


0.56 (3) 


0.57 (8) 


fb 


15.430 


0.69 (3) 


2.80 (5) 


residuals 




0.00134 





Table 5. Summary of the fixed and control parameters for the WD2004 
solution of the Kepler short-cadence light curve of KIC 10661783. For fur- 
ther details on the control pa rameters please see the WD2004 user guide 
IWilson & Van Hammejiooj) . 



Parameter WD2004name Star A StarB 



WD2004 operation mode 


MODE 




5 




Treatment of reflection 


MREF 


1 (simple) 


Number of reflections 


NREF 




1 




Third light 


el3 




0.0 




Effective temperature (K) 


TAVH, TAVC 


8000 




6500 


Rotation rates 


Fl,F2 


1.0 




1.0 


Gravity darkening 


GRl, GR2 


1.0 




0.3 


Bolometric albedos 


ALB 1 , ALB2 


1.0 




0.5 


Numerical accuracy 


Nl, N2 


60 




40 


Numerical accuracy 


NlL, N2L 


40 




30 


Bolometric linear LD coeff. 


XBOLl, XBOL2 


0.305 




0.110 


Bolometric nonlinear LD coeff. 


YBOLl, YBOL2 


0.416 




0.605 


Linear LD coefficient (starB) 


X2 






0.013 


Nonlinear LD coefficient 


Yl, Y2 


0.628 




0.717 



Table 6. Results of the WD2004 modelling process of the Kepler short- 
cadence light curve of KIC 10661783. The uncertainties come from the 
scatter of the individual solutions for the five subdivided datasets. 



Parameter WD2004 name Fitted value 



Star A potential 


PHSV 


2.625 ± 0.006 


Mass ratio 


RM 


0.0626 ± 0.0008 


Orbital inclination (°) 


XINCL 


82.80 ±0.18 


Bolometric albedo for star A 


ALB 1 


2.69 ±0.12 


Bolometric albedo for star B 


ALB 1 


1.09 ±0.06 


Light from star A 


HLUM 


0.13611 ±0.00004 


Light from star B 


CLUM 


0.01 168 ±0.00010 


Linear LD coefficient star A 


Xl 


0.108 ±0.006 


Fractional radius of star A 




0.3992 ±0.0011 


Fractional radius of star B 




0.1787 ±0.0007 



© 0000 RAS, MNRAS 000, 000-000 



KIC 10661 783: an eclipsing 5 Scuti star 9 



T 




5.10 - 



5.15 



-0.2 



0.0 



0.2 

Orbital phase 



0.4 



0.6 



Figure 6. Best WD2004 fit (red line) to the phased and binned Kepler short-cadence light curve of KIC 10661783 (blue points). The residuals of the fit are 
offset from zero to appear at the base of the plot. 



orbital motion. The gravity darkening coefficients were set to 1.0 
for the hotter star, which has a radiative atmosphere, and 0.3 for 
the co oler star, which probably has a convective atmosphere (Claret 
2000). Changes in these values do not have much effect on the light 
curve solution. 

The albedos were set to 1.0 and 0.5 respectively, as expected 
for stellar atmospheres of this type. The T c s s of the stars were put 
to 8000 K and 6500 K and were not adjusted during the fitting pro- 
cess. We instead fitted directly for the light contributions of the two 
stars. The T c s s are therefore used only in the specification the ini- 
tial limb darkening (LD) coefficients, and do not affect the shape 
of the calculated light curves. Th e square-root LD law was adopted 
jKlinglesmith &Sobiesldll97Cl) and initial coefficients were taken 
from the tabulations of IVan Hammd d 1993b . The nonlinear LD co- 
efficients were fixed and the line ar LD coefficient of star A was 
included as a parameter of the fit (ISouthwo rth 2008, 2010). Other 
fitted parameters were the potentials of the two stars, their mass 
ratio and the orbital inclination. 

Our initial attempts to fit the light curve were not very suc- 
cessful, which is unsurprising given the breathtaking quality of 
the Kepler light curve. Difficulties centred around matching the 
eclipses and the outside-eclipse variations simultaneously. At- 
tempts to solve this by fitting for rotational distortions or the grav- 
ity brightening exponents of the stars did not help significantly. We 
were finally able to obtain a much improved match to the data by 
including the stellar albedos as fitted parameters. The drawback to 
this approach is that the resulting value for the primary star is sub- 
stantially greater than unity, implying that the reflected/reprocessed 
light from the stellar surfaces is more than the incident light. We 
caution that this should not be taken as evidence for the sponta- 
neous creation of energy or the possibility of a perpetual motion 
machine. We attribute this phenomenon to the imperfection of the 
physics included in eclipsing binary light curve models. 

It may also be the case that the pulsating nature of star A is 
affecting the light curve more than we expect. The pulsations were 
removed from the light curve using their mean amplitudes, whereas 
the light contribution from star A is variable with orbital phase. The 



pulsation amplitudes will be more representative of the outside- 
eclipse brightness level, and will be overestimated during primary 
eclipse (when parts of the pulsating star are masked by its compan- 
ion) and underestimated during secondary eclipse (when the light 
from the pulsating star is no longer diluted by its colleague). A 
plot of the binned light curve versus the best WD2004 fit is shown 
in Fig.[6] Residual pulsation effects are certainly visible as 'wavi- 
ness' during the totality of secondary eclipse. They are probably 
also existe nt through primary ec lipse but are less obvious to the 
naked eye. Mac eroni et alj f20091 presented a study of HD 174884, 
a late-B-type EB observed by the CoRoT satellite. These authors 
found short-period correlations in the residuals of their fit to the 
light curve, and attributed these to pulsations which were phase- 
locked to 8 and 13 times the orbital frequency. Such a phenomenon 
may occur for KIC 10661783, and would result in pulsations which 
do not average down when the light curve is phase-binned. 

We finally arrived at two best-fitting light curve models, corre- 
sponding to two different binary configurations. The first is that of 
a detached binary system (WD2004 'mode 0'), and the other is for 
the case when the secondary star is filling its Roche lobe ('mode 
5'). The quality of the fit it essentially the same for both, so they 
are statistically indistinguishable with the current photometric data. 
The main difference between the two solutions is in the mass ratio, 
which is 0.25 for the former and 0.06 for the latter. Based on pre- 
liminary spectroscopic results we are able to reject the former solu- 
tion, which means that KIC 10661783 is a semi-detached eclipsing 
binary system. 

The final characteristics of the semi-detached eclipse model 
are given in Table[5](fixed and control parameters) and Table|6](fit- 
ted parameters). The best fit is shown in Fig. [6] for which the root- 
mean-square of the residuals is 0.95 mmag). Future work on this 
EB would benefit from a more extensive light curve, as this would 
allow the pulsational effects to be better averaged out as well as 
more accurately modelled during the eclipse phases. 

WD2004 returns error estimates calculated from the covari- 
ance matrix evaluated around the best fit to the light curve. These 
error estimates are known to become optimistic in the presence of 
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correlations between parameters (e.g. Maceroni & Rucinskill 19971 ; 
iBruntt et alj|2006t IPavlovski et al.|[2009fT In order to obtain more 
realistic error estimates we sequentially sliced up the short-cadence 
light curve into five different datasets, which were individually 
phase-binned and fitted with the WD2004 code. The error estimates 
were then calculated as the standard deviation of the five different 
values for each fitted parameter. We did not divide the standard de- 
viations by the square root of the number of estimates for each pa- 
rameter as we are in the regime where systematic errors dominate 
the fit (Fig.[6]l, and these are not averaged down by the provision 
of more data. Given that we were not able to obtain a perfect fit 
to the data, there is a systematic error on top of the quoted error- 
bars which is difficult to assess before improved light curve models 
become available. 



6 SUMMARY AND CONCLUSIONS 

We have presented Kepler photometry of KIC 10661783, a semi- 
detached binary system which shows total eclipses and multi- 
periodic <5 Scuti pulsations. The Kepler short-cadence and long- 
cadence data were augmented with SuperWASP observations in or- 
der to precisely measure the orbital period of the system. The short- 
cadence data comprise 38 000 datapoints obtained over 27 days and 
with a duty cycle of 89%. These were fitted with the JKTEBOP code 
in order to obtain a good morphological match to the light curve. 
This fit was then subtracted, leaving behind the pulsation signa- 
tures. 

This residual light curve was subjected to a frequency anal- 
ysis which revealed at least 68 frequencies of which 55 or more 
can be attributed to pulsation modes and the rest to harmonics of 
the orbital frequency. The frequency analysis was curtailed at this 
point due primarily to the limited frequency resolution of the data. 
The remaining unmodelled variation in the residual light curve is 
mostly pulsations whose frequencies are too closely spaced to be 
resolved. The main frequency range of the variability signal lies 
between 18 and 31 d _1 , and possesses amplitudes between 0.1 and 
4 mmag. One harmonic term (2 x /) and a few combination fre- 
quencies (/; + fj) have been detected. A dataset covering a longer 
time interval will enable the detection of hundreds of pulsation fre- 
quencies. From plots of the residuals of the 68-frequency fit to the 
data, we attribute the pulsation activity to the primary star. 

We then removed the pulsation signatures from the short- 
cadence light curve, leaving behind the variations due to binarity 
(eclipses, reflection effect and ellipsoidal effect). The light curve 
was phase-binned into 230 points, with a denser sampling for the 
eclipse phases than other phases, and modelled using the Wilson- 
Devinney code. We were able to get a reasonable but not com- 
pletely satisfactory fit to the binned data, and then only by push- 
ing the albedo of the primary star to a higher value than physically 
expected. We diagnose a problem with the approximate treatment 
of reprocessing and the reflection effect, which becomes apparent 
for light curves of the spectacular quality now attained by the Ke- 
pler satellite. Our fit to the binned light curve is also compromised 
by the residual pulsation signal apparent in the data. This arises be- 
cause the contribution of light from the pulsating primary star - and 
therefore the pulsation amplitudes - is variable throughout the orbit 
(most obviously during primary eclipse). 

After finding the best fit to the light curve we obtained error 
estimates for the fitted parameters by splitting the data into five sub- 
sets and modelling them individually. The fractional radii of the two 
stars (their radii expressed as a fraction of the semimajor axis) are 



measured to precisions of 0.3% and 0.4%, respectively. We have al- 
ready begun to obtain spectroscopic radial velocity observations of 
the two stars which, when combined with the photometric model, 
will yield direct measurements of the masses and radii of the two 
stars. The spectra will also be useful for measuring the T e gs and 
chemical abundances of our target. Although the secondary com- 
ponent contributes only 8.5% of the total light of the system, its 
atmospheric characteristics should be easily accessible using spec- 
tral disentangling techniques. We are already obtaining additional 
Kepler photometric observations of KIC 10661783, from which we 
will be able to measure a larger and more precise set of pulsation 
frequencies. These results could then be used to identify the pulsa- 
tion modes of the primary star. 

The semi-detached nature of KIC 10661783 means that it 
is a member of the class of oEA (oscillating Algol) stars 
dMkrtichian et al. 2002, 2003). The nu mber of confirmed oEA sys- 
tems is currently only about twenty dRodriguez et al.l |2010h. of 
which one is a high-amplitude S Scuti pulsator dChristiansen et al.l 
l2007h . The difficulty in performing asteroseismology from the 
ground means that the largest number of pulsation frequencies 
detected in an oEA system was previously only eight (YCam 
iRodriguez etai]|20ld) . KIC 10661783 is therefore the oEA with 
the richest pulsation spectrum, making it a very promising ob- 
ject for improving our understanding of the physical properties of 
5 Scuti stars. 
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